Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap.

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

sRNA-Effector: A tool to expedite discovery of small RNA regulators.

microRNAs (miRNAs) are small regulatory RNAs that repress target mRNA transcripts through base pairing. Although the mechanisms of miRNA production and function are clearly established, new insights into miRNA regulation or miRNA-mediated gene silencing are still emerging. In order to facilitate the discovery of miRNA regulators or effectors, we have developed sRNA-Effector, a machine learning algorithm trained on enhanced crosslinking and immunoprecipitation sequencing and RNA sequencing data following knockdown of specific genes. sRNA-Effector can accurately identify known miRNA biogenesis and effector proteins and identifies 9 putative regulators of miRNA function, including serine/threonine kinase STK33, splicing factor SFPQ, and proto-oncogene BMI1. We validated the role of STK33, SFPQ, and BMI1 in miRNA regulation, showing that sRNA-Effector is useful for identifying new players in small RNA biology. sRNA-Effector will be a web tool available for all researchers to identify potential miRNA regulators in any cell line of interest.

Integrative analysis of DNA replication origins and ORC/MCM binding sites in human cells reveals a lack of overlap.

Based on experimentally determined average inter-origin distances of ∼100 kb, DNA replication initiates from ∼50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the Origin Recognition Complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and 5 ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ∼7.5 million union origins identified by all datasets, only 0.27% were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques (20,250 shared origins), suggesting extensive variability in origin usage and identification. 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF binding sites, G-quadruplex sites and activating histone marks, these overlaps are comparable or less than that of known Transcription Start Sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ∼13,000 reproducible ORC binding sites in human cancer cells, and only 4.5% were within 1 kb of the ∼11,000 union MCM2-7 binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, S. cerevisiae . Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

XRN2 suppresses aberrant entry of tRNA trailers into argonaute in humans and Arabidopsis.

MicroRNAs (miRNAs) are a well-characterized class of small RNAs (sRNAs) that regulate gene expression post-transcriptionally. miRNAs function within a complex milieu of other sRNAs of similar size and abundance, with the best characterized being tRNA fragments or tRFs. The mechanism by which the RNA-induced silencing complex (RISC) selects for specific sRNAs over others is not entirely understood in human cells. Several highly expressed tRNA trailers (tRF-1s) are strikingly similar to microRNAs in length but are generally excluded from the microRNA effector pathway. This exclusion provides a paradigm for identifying mechanisms of RISC selectivity. Here, we show that 5' to 3' exoribonuclease XRN2 contributes to human RISC selectivity. Although highly abundant, tRF-1s are highly unstable and degraded by XRN2 which blocks tRF-1 accumulation in RISC. We also find that XRN mediated degradation of tRF-1s and subsequent exclusion from RISC is conserved in plants. Our findings reveal a conserved mechanism that prevents aberrant entry of a class of highly produced sRNAs into Ago2.

Depletion of the m1A writer TRMT6/TRMT61A reduces proliferation and resistance against cellular stress in bladder cancer.

Background: Bladder cancer (BLCA) is a common and deadly disease that results in a reduced quality of life for the patients and a significant economic burden on society. A better understanding of tumorigenesis is needed to improve clinical outcomes. Recent evidence places the RNA modification m1A and its regulatory proteins TRMT6/TRMT61A and ALKBH3 in BLCA pathogenesis. Methods: TRMT6/TRMT61A, ALKBH1, and ALKBH3 expression was examined in human BLCA cell lines and a normal urinary tract epithelium cell line through qRT-PCR and western blot analysis. Prestoblue Cell Viability Reagent, wound-healing assay, and live-cell imaging-based cell displacement analysis, were conducted to assess proliferation, migration, and displacement of this BLCA cell line panel. Cell survival was assessed after inducing cellular stress and activating the unfolded protein response (UPR) with tunicamycin. Moreover, siRNA-mediated gene silencing in two BLCA cell lines (5637 and HT1197) was conducted to investigate the biological roles of TRMT6/TRMT61A. Results: Heterogeneous morphology, proliferation, displacement, tunicamycin sensitivity, and expression levels of m1A regulators were observed among the panel of cell lines examined. In general, TRMT61A expression was increased in BLCA cell lines when compared to SV-HUC-1. Depletion of TRMT6/TRMT61A reduced proliferation capacity in both 5637 and HT1197 cell lines. The average cell displacement of 5637 was also reduced upon TRMT6/TRMT61A depletion. Interestingly, TRMT6/TRMT61A depletion decreased mRNA expression of targets associated with the ATF6-branch of the UPR in 5637 but not in HT1197. Moreover, cell survival after induction of cellular stress was compromised after TRMT6/TRMT61A knockdown in 5637 but not in HT1197 cells. Conclusion: The findings suggest that TRMT6/TRMT61A plays an oncogenic role in BLCA and is involved in desensitizing BLCA cells against cellular stress. Further investigation into the regulation of TRMT6/TRMT61A expression and its impact on cellular stress tolerance may provide insights for future BLCA treatment.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow.

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

Cryo-EM structure of an extracellular Geobacter OmcE cytochrome filament reveals tetrahaem packing.

Electrically conductive appendages from the anaerobic bacterium Geobacter sulfurreducens were first observed two decades ago, with genetic and biochemical data suggesting that conductive fibres were type IV pili. Recently, an extracellular conductive filament of G. sulfurreducens was found to contain polymerized c-type cytochrome OmcS subunits, not pilin subunits. Here we report that G. sulfurreducens also produces a second, thinner appendage comprised of cytochrome OmcE subunits and solve its structure using cryo-electron microscopy at ~4.3 Å resolution. Although OmcE and OmcS subunits have no overall sequence or structural similarities, upon polymerization both form filaments that share a conserved haem packing arrangement in which haems are coordinated by histidines in adjacent subunits. Unlike OmcS filaments, OmcE filaments are highly glycosylated. In extracellular fractions from G. sulfurreducens, we detected type IV pili comprising PilA-N and -C chains, along with abundant B-DNA. OmcE is the second cytochrome filament to be characterized using structural and biophysical methods. We propose that there is a broad class of conductive bacterial appendages with conserved haem packing (rather than sequence homology) that enable long-distance electron transport to chemicals or other microbial cells.

tRForest: a novel random forest-based algorithm for tRNA-derived fragment target prediction.

tRNA fragments (tRFs) are small RNAs comparable to the size and function of miRNAs. tRFs are generally Dicer independent, are found associated with Ago, and can repress expression of genes post-transcriptionally. Given that this expands the repertoire of small RNAs capable of post-transcriptional gene expression, it is important to predict tRF targets with confidence. Some attempts have been made to predict tRF targets, but are limited in the scope of tRF classes used in prediction or limited in feature selection. We hypothesized that established miRNA target prediction features applied to tRFs through a random forest machine learning algorithm will immensely improve tRF target prediction. Using this approach, we show significant improvements in tRF target prediction for all classes of tRFs and validate our predictions in two independent cell lines. Finally, Gene Ontology analysis suggests that among the tRFs conserved between mice and humans, the predicted targets are enriched significantly in neuronal function, and we show this specifically for tRF-3009a. These improvements to tRF target prediction further our understanding of tRF function broadly across species and provide avenues for testing novel roles for tRFs in biology. We have created a publicly available website for the targets of tRFs predicted by tRForest.

DeepTracer-ID: De novo protein identification from cryo-EM maps.

The recent revolution in cryo-electron microscopy (cryo-EM) has made it possible to determine macromolecular structures directly from cell extracts. However, identifying the correct protein from the cryo-EM map is still challenging and often needs additional sequence information from other techniques, such as tandem mass spectrometry and/or bioinformatics. Here, we present DeepTracer-ID, a server-based approach to identify the candidate protein in a user-provided organism de novo from a cryo-EM map, without the need for additional information. Our method first uses DeepTracer to generate a protein backbone model that best represents the cryo-EM map, and this model is then searched against the library of AlphaFold2 predictions for all proteins in the given organism. This method is highly accurate and robust for high-resolution cryo-EM maps: in all 13 experimental maps tested blindly, DeepTracer-ID identified the correct proteins as the top candidates. Eight of the maps were of known structures, while the other five unpublished maps were validated by prior protein annotation and careful inspection of the model refined into the map. The program also showed promising results for both homomeric and heteromeric protein complexes. This platform is possible because of the recent breakthroughs in large-scale three-dimensional protein structure prediction.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer.

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

Spindle-shaped archaeal viruses evolved from rod-shaped ancestors to package a larger genome.

Spindle- or lemon-shaped viruses infect archaea in diverse environments. Due to the highly pleomorphic nature of these virions, which can be found with cylindrical tails emanating from the spindle-shaped body, structural studies of these capsids have been challenging. We have determined the atomic structure of the capsid of Sulfolobus monocaudavirus 1, a virus that infects hosts living in nearly boiling acid. A highly hydrophobic protein, likely integrated into the host membrane before the virions assemble, forms 7 strands that slide past each other in both the tails and the spindle body. We observe the discrete steps that occur as the tail tubes expand, and these are due to highly conserved quasiequivalent interactions with neighboring subunits maintained despite significant diameter changes. Our results show how helical assemblies can vary their diameters, becoming nearly spherical to package a larger genome and suggest how all spindle-shaped viruses have evolved from archaeal rod-like viruses.

ATAC-Seq-based Identification of Extrachromosomal Circular DNA in Mammalian Cells and Its Validation Using Inverse PCR and FISH.

Recent studies from multiple labs including ours have demonstrated the importance of extrachromosomal circular DNA (eccDNA) from yeast to humans ( Shibata et al., 2012 ; Dillon et al., 2015 ; Møller et al., 2016 ; Kumar et al., 2017 ; Turner et al., 2017 ; Kim et al., 2020 ). More recently, it has been found that cancer cells obtain a selective advantage by amplifying oncogenes on eccDNA, which drives genomic instability ( Wu et al., 2019 ; Kim et al., 2020 ). Previously, we have purified circular DNA and enriched the population using rolling circle amplification followed by high-throughput sequencing for the identification of eccDNA based on the unique junctional sequence. However, eccDNA identification by rolling circle amplification is biased toward small circles. Here, we report a rolling circle-independent method to detect eccDNA in human cancer cells. We demonstrate a sensitive and robust step-by-step workflow for finding novel eccDNAs using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) combined with a Circle_finder bioinformatics algorithm to predict the eccDNAs, followed by its validation using two independent methods, inverse PCR and metaphase FISH (Fluorescence in situ Hybridization).

Structural analysis of cross α-helical nanotubes provides insight into the designability of filamentous peptide nanomaterials.

The exquisite structure-function correlations observed in filamentous protein assemblies provide a paradigm for the design of synthetic peptide-based nanomaterials. However, the plasticity of quaternary structure in sequence-space and the lability of helical symmetry present significant challenges to the de novo design and structural analysis of such filaments. Here, we describe a rational approach to design self-assembling peptide nanotubes based on controlling lateral interactions between protofilaments having an unusual cross-α supramolecular architecture. Near-atomic resolution cryo-EM structural analysis of seven designed nanotubes provides insight into the designability of interfaces within these synthetic peptide assemblies and identifies a non-native structural interaction based on a pair of arginine residues. This arginine clasp motif can robustly mediate cohesive interactions between protofilaments within the cross-α nanotubes. The structure of the resultant assemblies can be controlled through the sequence and length of the peptide subunits, which generates synthetic peptide filaments of similar dimensions to flagella and pili.

Structures of filamentous viruses infecting hyperthermophilic archaea explain DNA stabilization in extreme environments.

Living organisms expend metabolic energy to repair and maintain their genomes, while viruses protect their genetic material by completely passive means. We have used cryo-electron microscopy (cryo-EM) to solve the atomic structures of two filamentous double-stranded DNA viruses that infect archaeal hosts living in nearly boiling acid: Saccharolobus solfataricus rod-shaped virus 1 (SSRV1), at 2.8-Å resolution, and Sulfolobus islandicus filamentous virus (SIFV), at 4.0-Å resolution. The SIFV nucleocapsid is formed by a heterodimer of two homologous proteins and is membrane enveloped, while SSRV1 has a nucleocapsid formed by a homodimer and is not enveloped. In both, the capsid proteins wrap around the DNA and maintain it in an A-form. We suggest that the A-form is due to both a nonspecific desolvation of the DNA by the protein, and a specific coordination of the DNA phosphate groups by positively charged residues. We extend these observations by comparisons with four other archaeal filamentous viruses whose structures we have previously determined, and show that all 10 capsid proteins (from four heterodimers and two homodimers) have obvious structural homology while sequence similarity can be nonexistent. This arises from most capsid residues not being under any strong selective pressure. The inability to detect homology at the sequence level arises from the sampling of viruses in this part of the biosphere being extremely sparse. Comparative structural and genomic analyses suggest that nonenveloped archaeal viruses have evolved from enveloped viruses by shedding the membrane, indicating that this trait may be relatively easily lost during virus evolution.

Noncanonical Roles of tRNAs: tRNA Fragments and Beyond.

As one of the most abundant and conserved RNA species, transfer RNAs (tRNAs) are well known for their role in reading the codons on messenger RNAs and translating them into proteins. In this review, we discuss the noncanonical functions of tRNAs. These include tRNAs as precursors to novel small RNA molecules derived from tRNAs, also called tRNA-derived fragments, that are abundant across species and have diverse functions in different biological processes, including regulating protein translation, Argonaute-dependent gene silencing, and more. Furthermore, the role of tRNAs in biosynthesis and other regulatory pathways, including nutrient sensing, splicing, transcription, retroelement regulation, immune response, and apoptosis, is reviewed. Genome organization and sequence variation of tRNA genes are also discussed in light of their noncanonical functions. Lastly, we discuss the recent applications of tRNAs in genome editing and microbiome sequencing.

The structures of two archaeal type IV pili illuminate evolutionary relationships.

We have determined the cryo-electron microscopic (cryo-EM) structures of two archaeal type IV pili (T4P), from Pyrobaculum arsenaticum and Saccharolobus solfataricus, at 3.8 Å and 3.4 Å resolution, respectively. This triples the number of high resolution archaeal T4P structures, and allows us to pinpoint the evolutionary divergence of bacterial T4P, archaeal T4P and archaeal flagellar filaments. We suggest that extensive glycosylation previously observed in T4P of Sulfolobus islandicus is a response to an acidic environment, as at even higher temperatures in a neutral environment much less glycosylation is present for Pyrobaculum than for Sulfolobus and Saccharolobus pili. Consequently, the Pyrobaculum filaments do not display the remarkable stability of the Sulfolobus filaments in vitro. We identify the Saccharolobus and Pyrobaculum T4P as host receptors recognized by rudivirus SSRV1 and tristromavirus PFV2, respectively. Our results illuminate the evolutionary relationships among bacterial and archaeal T4P filaments and provide insights into archaeal virus-host interactions.

Structure of a filamentous virus uncovers familial ties within the archaeal virosphere.

Viruses infecting hyperthermophilic archaea represent one of the most enigmatic parts of the global virome, with viruses from different families showing no genomic relatedness to each other or to viruses of bacteria and eukaryotes. Tristromaviruses, which build enveloped filamentous virions and infect hyperthermophilic neutrophiles of the order Thermoproteales, represent one such enigmatic virus families. They do not share genes with viruses from other families and have been believed to represent an evolutionarily independent virus lineage. A cryo-electron microscopic reconstruction of the tristromavirus Pyrobaculum filamentous virus 2 at 3.4 Å resolution shows that the virion is constructed from two paralogous major capsid proteins (MCP) which transform the linear dsDNA genome of the virus into A-form by tightly wrapping around it. Unexpectedly, the two MCP are homologous to the capsid proteins of other filamentous archaeal viruses, uncovering a deep evolutionary relationship within the archaeal virosphere.

ATAC-seq identifies thousands of extrachromosomal circular DNA in cancer and cell lines.

Extrachromosomal circular DNAs (eccDNAs) are somatically mosaic and contribute to intercellular heterogeneity in normal and tumor cells. Because short eccDNAs are poorly chromatinized, we hypothesized that they are sequenced by tagmentation in ATAC-seq experiments without any enrichment of circular DNA. Indeed, ATAC-seq identified thousands of eccDNAs in cell lines that were validated by inverse PCR and by metaphase FISH. ATAC-seq in gliomas and glioblastomas identify hundreds of eccDNAs, including one containing the well-known EGFR gene amplicon from chr7. More than 18,000 eccDNAs, many carrying known cancer driver genes, are identified in a pan-cancer analysis of ATAC-seq libraries from 23 tumor types. Somatically mosaic eccDNAs are identified by ATAC-seq even before amplification is recognized by genome-wide copy number variation measurements. Thus, ATAC-seq is a sensitive method to detect eccDNA present in a tumor at the pre-amplification stage and can be used to predict resistance to therapy.

miR-206 family is important for mitochondrial and muscle function, but not essential for myogenesis in vitro.

miR-206, miR-1a-1, and miR-1a-2 are induced during differentiation of skeletal myoblasts and promote myogenesis in vitro. miR-206 is required for skeletal muscle regeneration in vivo. Although this miRNA family is hypothesized to play an essential role in differentiation, a triple knock-out (tKO) of the three genes has not been done to test this hypothesis. We report that tKO C2C12 myoblasts generated using CRISPR/Cas9 method differentiate despite the expected derepression of the miRNA targets. Surprisingly, their mitochondrial function is diminished. tKO mice demonstrate partial embryonic lethality, most likely due to the role of miR-1a in cardiac muscle differentiation. Two tKO mice survive and grow normally to adulthood with smaller myofiber diameter, diminished physical performance, and an increase in PAX7 positive satellite cells. Thus, unlike other miRNAs important in other differentiation pathways, the miR-206 family is not absolutely essential for myogenesis and is instead a modulator of optimal differentiation of skeletal myoblasts.